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We demonstrate a novel technique for the measurement of the charge carried by a colloidal particle. 
The technique uses the phenomenon of the resonance of a particle held in an optical tweezers trap 
and driven by a sinusoidal electric field. The trapped particle forms a strongly damped harmonic 
oscillator whose fluctuations are a function of 7, the ratio of the root-mean square average of the 
electric and thermal forces on the particle. At low applied fields, where 7 < 1, the particle is 
confined to the optical axis while at high fields (7 > 1) the probability distribution of the particle 
is double-peaked. The periodically-modulated thermal fluctuations are measured with nanometer 
sensitivity using an interferometric position detector. Charges, as low as a few elementary charges, 
can be measured with an uncertainty of about 0.25 e. This is significantly better than previous 
techniques and opens up new possibilities for the study of nonpolar suspensions. 



I. INTRODUCTION 

Although charge is often thought to play no role in 
nonpolar media there is ample [l|, 0, HI], if sometimes 
contradictory Q, evidence of surface charging in non- 
polar colloidal suspensions. Much of the reason for this 
uncertainty is the extremely low level of charge in non- 
polar environments. While a particle of radius a ~ 500 
nm might carry a surface charge eZ e R of the order of 
10 3 electrons in an aqueous dispersion (e r ~ 80), the 
charge in a nonpolar solvent will be some two or three 
orders of magnitude smaller. Although small in absolute 
terms, this charge still produces surprisingly strong ef- 
fects in low dielectric environments. For instance, using 
Coulomb's law the contact value of the interaction po- 
tential (in units of fcgT) between two colloidal spheres 
with radius a and charge eZ e g in a solvent of dielectric 
ratio e r is simply, 

k B T k B T 87re e r a \ 2a J ofr ' 1 ' 

where we have introduced the Bjerrum length A^ = 
e 2 / 1 ^TTeotrksT as a characteristic length scale of electro- 
static interactions. For dodecane, As is 28.0 nm so the 
repulsive interactions between 500 nm particles, carrying 
a typical charge of 10 electrons, is w 3 fcsT - big enough 
to have a dramatic effect on the structure of a colloidal 
suspension. 

The most direct route to monitor the extent of parti- 
cle charging is to measure the uniform drift velocity up 
in the presence of a weak electric field E. For a spher- 
ical particle, urj = fiE, where fi is the electrophorctic 
mobility. For an isolated colloidal sphere, in the low-salt 
regime, moving in a solvent of viscosity rj, Stokes' law im- 
plies fj, = Z e ge/ (6irr]a) since the ionic atmosphere around 
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the colloid is very diffuse. The low charges found in non- 
polar media therefore translate to electrophoretic mobili- 
ties some two or three orders of magnitude or so less than 
in water, at the same electric field. Rapid and repro- 
ducible measurements of such low mobilities presents a 
major challenge to conventional clcctrokinetic techniques 

a. 

Motivated by the growing acknowledgement of the im- 
portance of charge in nonpolar suspensions [j| and its 
technological significance [6[ , several novel methods have 
been developed to determine particle mobility in nonpo- 
lar suspensions. Schatzel and coworkers pioneered an 
innovative amplitude- weighted phase analysis signal pro- 
cessing scheme for the analysis of classical laser Dopplcr 
electrophoresis data. The technique, phase analysis light 
scattering (PALS), is significantly more sensitive than 
conventional laser Doppler electrophoresis. In practice 
however, collective particle motion induced by thermal 
convection or sedimentation restricts the achievable sen- 
sitivity of commercial PALS instruments to ~ 10 -10 m 2 
s _1 V -1 , which corresponds to colloid charges of ~ 10e 
on a particle of radius a — 500 nm. Video microscopy 
has been used by several authors. Perez and Lemaire [8j 
used a vibrating piezoelectric stage to null the oscillatory 
microscopic motion produced by an applied alternating 
field, while Strubbe et al. § have measured the diffu- 
sion constant and electrophoretic mobility of individual 
particles in a nonpolar suspension. Despite these devel- 
opments there is still no technique which can provide fast 
and accurate measurement of the extremely small levels 
of particle charge found in nonpolar suspensions. As a re- 
sult the microscopic charging mechanisms in non-aqueous 
media has remained largely problematic Q. Improve- 
ments in quantitative characterization are an essential 
first step to developing a more detailed understanding of 
the role of charge in nonpolar suspensions. 

In this paper, we present a simple technique that al- 
lows a direct measurement of the effective charge of col- 
loidal particles in nonpolar solvents. Our approach is 
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very different from standard clcctrophorctic light scatter- 
ing methods and uses the phenomenon of the driven reso- 
nance of a strongly-damped harmonic oscillator. We use 
the radiation pressure from a laser beam focused to a near 
diffraction-limit spot (an 'optical tweezer') to trap single 
colloidal particles in a dilute suspension. The focused 
beam produces a harmonic potential for the trapped par- 
ticle whose strength is proportional to the laser power. 
The hindered diffusion of the trapped colloid is modified 
by an oscillatory electric field. The modulated thermal 
fluctuations are measured with nanometer sensitivity and 
the charge extracted from the spectral density of the fluc- 
tuations at the field frequency. 

Our approach is similar in principle to the classical Mil- 
likan oil drop experiments [10 ] but uses light rather than 
gravity. In his famous experiment to prove the quantiza- 
tion of charge, Millikan measured the speed of a charged 
oil drop under gravity and its speed falling in a elec- 
tric field E. He switched on electric forces Fe = ZeE 
of between 300-1300 fN (1 fN = 10~ 15 N) on oil drops 
with radii between 2 and 6 /xm (each drop containing 
Z elementary charges, e) to counterbalance the gravita- 
tional force mg acting on the drop. Our experiments 
rely, instead, on using the random thermal forces leading 
to Brownian motion, to oppose an applied electric force. 
The balance of forces determines the diffusive motion of 
the particle, rather than the velocity of free-fall as in 
Mullikan's experiment. Below we demonstrate that, by 
using a weak optical trap and an interfcrometric position 
detection system, we can measure accurately a mean par- 
ticle charge of -2.9 e which, for the typical field strengths 
used, corresponds to an electric force of about 40 fN. We 
estimate that the sensitivity of our technique is about 
an order of magnitude smaller than this limit. We be- 
lieve our method is a significant improvement on existing 
techniques. It provides rapid, accurate and highly sen- 
sitive information on the status of particle charging in 
nonpolar suspensions. Related but less sensitive exper- 
iments have been described on highly-charged aqueous 
colloids by Garbow et al. [llj where the authors used a 
weakly-focused laser beam. This had the effect of push- 
ing particles along the optical axis, rather than stably 
confining them, as in our approach. 

This paper is organized as follows: After a description 
of the apparatus and measurement procedures, we give a 
short introduction to the theory of a driven overdamped 
harmonic oscillator before results obtained on a variety 
of different nonpolar suspensions are presented and dis- 
cussed in the final section. 



II. MATERIALS AND METHODS 

A. Materials. 

We use a model non-aqueous colloidal system of 
sterically-stabilized poly (methyl mcthacrylatc) spheres 
of radius a = 610 ± 30 nm suspended in dodecane. 



The particles were synthesized by dispersion polymer- 
ization, following the procedures detailed by Antl et al. 
[l2l ]. Electron microscopy revealed the particles were 
highly uniform in size with a radius polydispersity (root 
mean square variation / mean radius) of a = 0.046 
± 0.01. Dodecane was dried with activated molecular 
sieves (Acros, size 4A) and stored under dry nitrogen. 
Sodium bis(2-ethylhexyl) sulfosuccinate (NaAOT, Fluka 
BioChcmika Ultra 99%) was purified by dissolution in 
methanol and tumbled with activated charcoal. Zirconyl 
2-ethyl hexanoate [Zr(Oct)2] was purchased from Alfa 
Aesar (Heysham, UK) and used as received. 



B. Optical Tweezers. 

Figure Q] shows a schematic of our optical tweezer 
system. A fibre-coupled laser beam with wavelength 
A = 1064 nm (Ytterbium fiber laser, IPG Photonics, Ger- 
many) was focusscd to a diffraction-limited beam waist 
by a microscope objective (Plan-Neofluar, xlOO, N.A. 
1.3, Zeiss) mounted in an inverted microscope (Axiovert 
S100, Zeiss). The laser intensity (~ 10 mW) was var- 
ied using a A/2 waveplate and a polarizing beam-splitter 
cube placed in the beam path. The fluctuations in posi- 
tion of the trapped particle, with respect to the center of 
the trap, were measured with a quadrant photodetector 
(model QD50-4X, Centronics, UK) using back-plane op- 
tical interferometry [l3| . The particle positions were ac- 
quired by a Labview programme (National Instruments, 
Austin, Texas, USA) and digitized using a high-speed 
data acquisition card (National Instruments model PCI- 
MI0-16E-1) at 10 kHz. For the current experiments, the 
number of particle positions collected in each time trace 
was set at 2 18 , so that the duration of each measure- 
ment was « 26 s. Positions measured in voltages were 
converted into displacements in nanometers by recording 
the time-dependent mean-square displacement (Ai 2 (t)) 
of five particles from the same batch of particles, imme- 
diately before each set of electrophoresis measurements. 
The laboratory temperature was stabilized at 22 ± 1 °C. 
At this temperature the viscosity of dodecane is 1.378 x 
10 -3 Pa s and the dielectric constant is 2.002. 



C. Single Particle Optical Microelectrophoresis 
(SPOM). 

The microelectrophoresis cell was constructed from 
two platinum foil electrodes, 127 /im high and with a 
width of 2 mm, mounted in a cylindrical glass sample 
chamber, 1.13 mm high with a volume of ~ 90 [A. The 
cell was constructed in three stages. First, two platinum 
foil electrodes were carefully positioned parallel to each 
other onto a 5.0 cm x 2.5 cm rectangular glass cover- 
slip (170 fim thick), covered with a thick layer of a UV- 
activated glue (Loctite 350), using a low-magnification 
stereo microscope. The electrodes were securely fixed in 
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place by irradiation with a 100 W UV-lamp. Then a 1 
mm thick glass slide with an 1 cm diameter inset circular 
hole was placed on top of the electrodes and sealed to 
the electrodes with UV glue. Finally, immediately prior 
to use, the top of the chamber was closed by a 22 mm 
diameter 170 /im-thick circular coverslip. The use of op- 
tical quality cover-slips on both the bottom and top of 
the sample chamber minimized spherical abberation. 

The distance between the two electrodes was mea- 
sured, after assembly of the cell, using a calibrated micro- 
scope graticule as 128 /im ± The small gap between 
the electrodes makes it possible to generate electric fields 
on the order of tens of kilovolts per meter with just a 
few volts applied to the cell, which minimizes thermal 
instabilities due to Joule heating and avoids the electro- 
hydrodynamic fluid instabilities which occur at higher 
voltages [l3 |. The electric field was estimated from the 
expression E = XV/ d, where V is the applied voltage, d 
is the plate separation and, A is a factor correcting for 
the finite size of the electrodes. Analytical calculations 
[l5T ]. for the field on the midplane equidistant between 
two semi-infinite plates, give A = 0.969 for the rectangu- 
lar electrodes used. Accumulation of ions near the elec- 
trode and the generation of a non-uniform space charge 
was suppressed by using an ac voltage. Electrosmosis was 
minimized by positioning the electrodes 60 /jm above the 
bottom glass wall, so that the tangential electric field at 
the wall was reduced. Moving the position of the trapped 
particle away from the electrode center by 10 ^m, in any 
direction, changed the recorded charge by less than 2.5%, 
confirming the absence of electrosmotic fluid flows. 

The sample chamber was filled with a dilute suspension 
of colloidal particles with a volume fraction of ~ 10~ 5 . 
The cell was carefully centered so that the laser passed 
through a plane equidistant between the two electrodes 
and the focus was midway between the top and bottom 
surfaces. The particle was held at least 65 /im from 
the nearest surface, to ensure Stokes' law was applicable. 
The charge on the particle was analysed using a purpose- 
written software package written in IDL (Research Sys- 
tems, Boulder, Colorado). The sign was determined by 
reducing the frequency of the applied electric field to 0.5 
Hz, blocking the laser beam so that the particle was mo- 
mentarily released from the optical trap, and following 
the oscillatory motion of the free particle on a CCD cam- 
era. 

The accuracy of SPOM was checked by measuring 
the mobility of a 850 nm (radius) PMMA particle in 
dodecane (with 100 mM NaAOT) using a commercial 
PALS instrument (Brookhaven ZetaPlus), ^pals, and 
the SPOM technique, /xspom- The good agreement be- 
tween //spom ((-5.4± 0.5) x 10" 10 m 2 s" 1 V" 1 ) and 
Mpals ((-6.0± 0.5) x 10~ 10 m 2 s" 1 V" 1 ) confirms that 
the optical microelectrophoresis technique yields accu- 
rate particle mobility data. 



III. THEORY 

A colloidal particle illuminated by a laser beam en- 
counters a gradient force that is oriented in the direction 
of the intensity gradient and a scattering force that is ori- 
ented in the direction of the incident light. In a tightly- 
focused beam the gradient forces dominate and lead to 
a strong restoring force, which causes the particle to be 
confined near the focus of the beam. To a first approx- 
imation, the particle can be modeled as a mass m in a 
three-dimensional harmonic potential. The potential is 
characterized by two spring constants, one in the axial 
and one in the radial direction. Here, for simplicity, we 
analyze only the one-dimensional motion along the radial 
coordinate x, where the corresponding force constant is 
k H . 

Placed in an external electric field, the total force on 
a charged particle is a sum of three contributions: a har- 
monic force, —kHx(t), arising from the optical trap; a 
random force 7Z(t) that represents thermal forces at the 
temperature T; and an external periodic force F p (t) = 
Asm(uj p t + 4>) which is characterized by the frequency 
Up, amplitude A and initial phase 4> of the external field. 
Defining the effective charge on the particle as Z c g (in 
units of the fundamental charge e) then the amplitude of 
the oscillatory force is 

A = Z cS eE, (2) 

where E is the maximum electric field strength at 
the particle center. In a single particle optical mi- 
croelectrophoresis experiment, the quantity ultimately 
measured is the time-dependent particle autocorrelation 
function C (r) , or its Fourier-transform, the spectral den- 
sity Below we derive expressions for both quanti- 
ties, in terms of the effective charge Z e ff. 

A. A periodically-driven Brownian Oscillator 

The fluctuations of a periodically-driven Brownian 
oscillator are described, in the conventional Ornstein- 
Uhlenbeck theory by the Langevin equation 

mx(t) + £x(t) + k H x{t) = Asm(uj p t + <t>) + TZ{t). (3) 

Here x{t) is the particle trajectory and £ = Gnija is the 
friction coefficient of a sphere of radius a moving in a 
medium of viscosity 77. The random thermal forces are 
modeled by a Gaussian process lZ(t) with the moments, 

(K(t)) = 0, and {K(t)K(t')) = 2£k B T5{t - t'). 

. (4) 

This definition ensures that the Brownian particle is in 
thermal equilibrium at a temperature T, in the absence 
of any driving field i.e. for A = 0. 

A colloidal particle, once moving, loses momentum 
rapidly as a result of viscous losses. The characteris- 
tic time for this decay is ts = rn/S,, which is about 
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10 3 times smaller than the time resolution of our experi- 
ments, when digitizing at a 10kHz sampling rate. Conse- 
quently the inertial terms may safely be dropped so that 
Eq. [3] now reads 

£x(t) + k H x(t) = Asm(u p t + (/))+ K(t). (5) 

The equation is linear so that the general solution may 
be written fl7l. 



x(t) = x n (t) + x p (t), 



(6) 



where x n is the solution in the presence of random ther- 
mal forces only (i.e. A = 0) and x p is the solution when 
only periodic forces act (i.e. T = 0). We now consider 
each solution in turn. 

The driven oscillator is deterministic in the absence of 
any random thermal noise. The solution is well known 



any 



x p (t) = 



A 



k H [1 + (ujp/uj c ) 



— sm(ujpt + (j)- A). (7) 



2lV 



The particle motion lags the driving field by a constant 
phase factor, 



A = tan 



-l "-p 



(8) 



which varies with the relative sizes of the driving and the 
corner frequency, 



w c = k H /£. 



(9) 



The frequency ratio lo p /uj c also controls the amplitude of 
the particle motion. For low frequencies (uj p <C w c )j the 
amplitude is essentially frequency independent. While, 
at high frequencies (oj p ^> ui c ) the response decreases 
rapidly with increasing u> p . In between, Eq. [7] reveals 
a monotonic oscillation spectrum, without the dynamic 
resonance at w p = uj c , characteristic of underdamped sys- 
tems. The relative phase of the particle motion depends 
on the initial phase <f> of the applied force. In our ex- 
periments <f> is a fluctuating stochastic variable with a 
coherence time t p , which depends on the stability of the 
voltage generator. The phase is essentially randomized 
on times t 3> t p so time and phase averages are equiv- 
alent. Multiplying the trajectories at t and t' and aver- 
aging over 0, assuming the initial value of <fi is uniformly 
distributed over the interval [0, 2tt], yields the correlation 
function, 



(xp(t)x p (t')), = 



A 2 



■coaujp{t - t'), (10) 



* 2k\ [1 + (cj p /u c y 

where the subscript denotes the phase-averaging. Since 
the motion is totally deterministic, (x p (t)x p (t')} ^ oscil- 
lates continuously between positive and negative values. 

The dynamics of the trapped particle is a linear super- 
position of periodic and random motion. The correlation 



function of a pur ely Brownian oscillator is given by Doi 
and Edwards [l8| as, 



(x n (T)x n (0)) 



k R T 



exp(-w c r). 



(11) 



The mean-square displacement is accordingly (xfy^ — 
ksT/kH and, since F = —knx, the mean-square thermal 
forces on the trapped particle are (i 7 ^)^ = ksTkn- 

B. Correlation function and spectral density 

The correlation function c(r) of the periodically-driven 
Brownian oscillator is defined by 



C (T = t'-t) = ((x(t)x(t')) R ),, 



(12) 



where the averages arc taken over both noise and phase. 
Since the thermal and periodic forces are uncorrelated, 
c(t) is simply a sum of the correlation functions for pe- 
riodic and purely thermal motion, 



c r 



k B T 



exp(— lj c t) + 



,4 2 



2k% [1 4- {u}p/uj c f 



(13) 

It is more convenient to work with the normalized func- 
tion C(t) = c(t)/c(0), which may be written as 



C(r) 



1 



i + r 



■ exp(-u) c r) 



T 



14-7 



2 cosw p t, (14) 



where we have introduced 7 2 , the scaled ratio of the 
mean-square periodic and Brownian forces, 



7 



1 4- {U!p/Uj c ) 2 



(15) 



Here <F p 2 ) = Z 2 H e 2 E 2 /2 and (F 2 ) = k B Tk H . The force 
ratio 7 appears naturally in the theory of a driven Brow- 
nian oscillator and, as we show below, is also the quantity 
most readily extracted from experiment. Knowing 7, the 
charge on the particle follows immediately as, 



Vt?^^- (16) 



This expression simplifies at low frequencies (lu p <C to c ) 
to, 



e\Z eS \ = —y / 2k B Tk H 
hi 



(17) 



which does not depend on the radius of the trapped par- 
ticle. 

Figure [3] shows the force ratio 7 calculated for typical 
experimental parameters. In the weak-field limit, where 
7 <C 1, the motion of the trapped particle is dominated 
by random thermal forces. By contrast, in the strong- 
field limit where 7 3> 1, the thermal forces arc only a rel- 
atively small perturbation and the electrophoretic forces 
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dominate. The cross-over between these two regimes oc- 
curs at a particle charge of Z e g ~ 20, for typical field 
strengths. At this point the electric forces on the parti- 
cle are a few hundred femtonewtons. 

The force ratio 7 may be determined experimentally 
from the correlation function C(r) or, equivalently, from 
the spectral density, 

1 Z 100 

I(Q,) = — c(r) exp(-iOr)dr. (18) 

As we discuss below there are advantages to both ap- 
proaches. Inserting Eq. [13] into [18] reveals that the spec- 
tral density is a sum of two independent contributions, 

m) = ^ ^ + k -^fm - u p) + s{ n + . p)] , 

(19) 

a Lorentzian spectrum typical of Brownian motion in a 
harmonic potential, together with additional <5-spikes at 
the positive and negative field frequencies. Integrating 
the spikes over the frequency axis gives an expression for 
the mean-squared periodic displacement of the particle 
p 

/OO 7 rp 

/ (n)dn = -^ 7 2 . (20) 
-oo k H 

from which 7 can be found. 



C. Asymptotic probability distribution function 

We now turn to the probability distribution of the con- 
fined particle. Introducing the dimcnsionlcss position co- 
ordinate x = xj '(ksT '/fcf/) 1 / 2 , we define the probability 
of locating the Brownian particle in the region between 
x and x + dx as p(x)dx, where p(x) is the probability 
distribution function. 

Without a periodic force, the Brownian particle is at 
thermal equilibrium and x(t) follows the Boltzmann's dis- 
tribution, 

p n (x) = _l=exp(-a; 2 /2). (21) 

V Z7T 

Applying a time-dependent force perturbs the system 
and shifts the average distribution away from the equilib- 
rium result. When all transients have died out, the prob- 
ability approaches a new asymptotic distribution which 
is a periodic function of the time t. As discussed ear- 
lier, the only experimentally accessible quantities are the 
phase-averaged functions. Analytical expressions for the 
phase-avera ged distributions have been proposed by Jung 
and Hanggi |l9j, for arbitrary friction coefficients. In the 
overdamped limit, their result reduces to the infinite se- 
ries, 



p(x) 



exp 



x 2 + 7 2 " 



Io (\/2yc) + 2^(-l) fe 4. f y J hk (V2 7 x) 



k=l 



(22) 



r 



where Ik (x) is a modified Bessel function [2(j ■ Figure d] 
shows the effect of a finite external field on the probability 
distribution. With increasing amplitude, the probability 
p(x) first broadens until, at 7 ~ 7*, the distribution be- 
comes bimodal and develops two peaks around x ~ ± 7. 
Further increases in the external field lead the two peaks 
to split further apart, leaving a relatively flat distribution 
around the origin, as the particle becomes increasingly 
localized around the two extremities. 

The changes in the probability distribution with exter- 
nal field can be clearly seen from the curvature of p(x) 
near to the origin. Expanding Eq. 1221 around x = gives, 
to terms in x 2 , 



p{x) 



1 



exp 



7_ 
" 2 



(1 - kx 2 ) (23) 



where the curvature is 



k = — 
2 



1 - 7 2 + 7 2 



^i(7 2 /2) 



M7 2 /2) 



(24) 



k is positive for a single monomodal distribution and neg- 
ative for a bimodal distribution. The inset plot in Fig.H 
displays the variation of n with external field. The cur- 
vature first decreases rapidly with increasing 7, reflect- 
ing the broadening of p{x) evident at low modulation 
strengths, before changing sign for 7 > 7* = 1.257 as the 
distribution becomes bimodal. The small negative values 
of the curvature found at high fields correspond to the 
extremely flat shape of the distribution seen close to the 
origin at high 7. The curvature is a function solely of 
7 2 so that the distribution sharpens continuously with 
increasing frequency and does not show the resonant bi- 
modality evident in the underdamped limit [l9| . 

The width of the particle distribution is (a; 2 ), which 
from Eq. [T31 equals 



(25) 



where (■ ■ • ) denotes an average over both 1Z and <p. In 
the weak- field limit where 7 « 1. the width of the dis- 
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tribution approaches the thermal limit. ksT/kn- Broad- 
ening of the probability distribution is only significant 
in the strong- field limit (7 3> 1) where the width varies 
quadratically with the strength of the applied field. 



IV. RESULTS AND DISCUSSION 

A. Field and Frequency Dependence of Effective 
Charge 

To validate our technique, experiments were first car- 
ried out on a number of weakly-charged particles, as a 
function of the applied field strength and its frequency. 
The particles were taken from a dilute dispersion of 610 
nm (radius) poly (methyl mcthacrylate) particles in a 
0.035 wt% solution of a PHSA-g-PMMA copolymer in 
dodecane. Figure [5] shows a short portion of the x- and 
y-position of a trapped particle when a sinusoidal field is 
applied. The field was applied along the x-axis so peri- 
odic oscillations are expected in the x-signal and should 
be absent from the y-signal. Looking at the trajecto- 
ries plotted in Fig. El it is very difficult to see any sig- 
nificant differences between the x- and y-signals. The 
small clectrophoretic oscillations expected are 'buried' 
beneath much larger Brownian fluctuations and 7 < 1. 
To extract the small periodic excursions from the pre- 
dominantly random Brownian motion we determine the 
spectral distribution of the fluctuations. 

Figure [5] shows the normalized autocorrelation func- 
tions, C (r) , measured on the same particle, as the applied 
field was varied. The oscillatory component is clearly vis- 
ible, particularly at long times. Data collection is fast 
with each of the functions shown representing just 26 s 
worth of data. The alternating electric field had a fre- 
quency of 67 Hz and voltages of up to 10 V were ap- 
plied to the electrodes. The laser trap had a stiffness of 
kn = 4.7±0.f fN nm -1 and a corner frequency of 47 ± 1 
Hz. The agreement between the experimental data and 
the least-squares fits is excellent. Fig.[6£a) reveals that at 
short-times the correlation functions all essentially con- 
verge, even though their long time dependence is clearly 
very different. This is because on timescales shorter then 
the period of the applied field, the motion of the trapped 
particle is dominated by Brownian diffusion. Thermal 
fluctuations are independent of the field and, since they 
are uncorrelated, result in the exponential relaxation ap- 
parent at short times in C(r) (Eg. ITT]) . The oscillations 
evident at long times (Fig. [51(b)) reflect the periodic mo- 
tion produced by the applied field and accordingly grow 
as the applied field is increased. 

Fitting the measured correlation functions to the the- 
oretical expression for C(r) (Eq.[TJ} yields 7, the ratio of 
the rms electrophoretic and Brownian forces. The ratio 
is plotted in Figure [7] as a function of the applied field 
E. As expected, the experimental values of 7 are lin- 
early proportional to E, confirming that the charge on 
the particle is constant. The gradient of the plot yields 



an effective charge of \Z e g\ = 14.5 ± 0.2 e, equivalent to 
an electrophoretic mobility of /x = 1.46 ± 0.02 x 10~ 10 
m 2 s _1 V -1 . The charge on this particle is sufficiently 
low that, even at the maximum field strengths achiev- 
able, the system remains in the low- field limit (7 < 1) 
where Brownian forces dominate. 

For ultra-low charges where 7 <C 1 it becomes increas- 
ingly difficult to determine the oscillatory element in the 
autocorrelation function, because of fluctuations in the 
baseline. This can be seen in the data at 13.3 kV m _1 in 
Figure [6fb) where the oscillations are indistinct. To re- 
solve these low particle charges we perform a fast Fourier 
transform of the data and calculate the spectral den- 
sity. Figure [5] shows the field dependence of the resulting 
power spectra. The measured spectra are particularly 
simple, being a superposition of a Lorentzian spectrum 
(shown dashed) reflecting diffusive motion and a single 
(5-peak at the fundamental electrode drive frequency u> p , 
whose height increases with the applied field. The data 
contains no higher harmonics of u p confirming the lin- 
ear response of particle and field. To extract the power 
in the periodic signal, the data around the <5-spike was 
masked and the diffusive spectrum fitted to the first term 
of Eq. If 91 by adjusting the unknown corner frequency 
lu c . Subtracting the fitted Lorentzian from the measured 
power spectra yielded the signal spectrum. The mean- 
square electrophoretic displacement P S i g was obtained by 
integrating the signal spectrum around the peak at u> p . 
The resulting values for P S i g are plotted as a function of 
the square of the applied field, in the inset of Figure [8] 
As expected from Eq. [201 the mean-square displacement 
P S i g varies linearly with E 2 . The high positional sensi- 
tivity achievable is highlighted by the arrowed (5-peak in 
Fig. [8] which, while it remains clearly visible, corresponds 
to a mean-square displacement of only f .4 nm 2 . 

To explore how the frequency of the applied field in- 
fluenced the measurements we collected autocorrelations 
at a number of different field frequencies, taking care to 
ensure that the same particle was used in each case. The 
corner frequency of the trapped particle was uj c = 291 
rad s^ 1 and data was recorded at frequencies between 
25 rad s _1 and f930 rad s _1 , approximately an order of 
magnitude below and above the corner frequency of the 
optical trap. The amplitude of the applied field was fixed 
at 75.8 kV m -1 . The autocorrelations were analyzed as 
described above to determine 7 2 at each frequency. The 
results are summarized in Figure [9] As expected, from 
the discussion in 1111 A[ the response of the trapped par- 
ticle to the applied field decreases with increasing fre- 
quency. This drop in sensitivity is particularly sharp 
when the driving frequency exceeds the corner frequency 
of the optical trap. Plotting 7 2 as a function of the fre- 
quency ratio 1 + (u>p/u! c ) 2 confirms that 7 2 is inversely 
dependent on 1 + (w p /cj c ) 2 , in agreement with the predic- 
tions of Eq. 1151 The highly linear dependency evident in 
Fig. O confirms the validity of the driven harmonic model 
and indicates that the charge on the particle is frequency 
independent. Analysis gives the surface charge on the 



7 



particle as \Z eS \ = 15.8 ± 0.2. 



Field Dependence of Probability Distribution 
Function 



So far, we have found excellent agreement between 
our measurements and the theoretical predictions for a 
periodically-driven oscillator. The comparison has how- 
ever been limited to weakly-charged particles where 7 < 
1. To extend our understanding, we now look at more 
highly-charged particles where strong electrophoretic 
forces dominate the weaker Brownian forces. Hsu et al. 
has shown Q that non-polar suspensions may be charged 
by adding the surfactant, sodium bis(2-ethylhexyl) sul- 
fosuccinate (Na-AOT) , which forms reverse micelles in a 
non-polar solvent. Adding Na-AOT to our suspensions 
at concentrations of 100 mM increases the average charge 
(in electrons) on the PMMA particles by a factor of four 
to Z e g w —50. Using these more highly-charged parti- 
cles we are able to investigate the high-field limit where 
7 > 1. 

Figure [TD] shows the field dependence of the parti- 
cle probability. The normalized distribution function 
is plotted as a function of the dimensionlcss position 
x = xjiksT I^h) 1 ^ 2 i where the re-location of the sphere 
has been scaled by the size of the corresponding rms 
Brownian fluctuations in the optical trap. To correct for 
small asymmetries in the trap caused by optical misalign- 
ment, the measured probability was averaged together at 
each positive and negative displacement. The resulting 
symmetrized distributions, p s (x) = \p(x) + p(—x)]/2, are 
plotted in Figure [TOl At zero field the probability distri- 
bution is a simple Gaussian function, as expected from 
Eq. [5T] Increasing the field to E « 23 kV m _1 causes 
the distribution to broaden, although it remains singly- 
peaked. At fields of 46 kV m _1 and above there is a 
dramatic change in p s (x). The central peak in the dis- 
tribution first splits into two broad maxima, before then 
moving further apart in x as the field is increased still 
further. Using the experimentally-determined charge on 
the particle (Z c g = — 48.3 ± 0.3) and the measured opti- 
cal trap stiffness we calculate the force ratio 7 from the 
applied electric field E. In this way we have confirmed 
that the probability distribution first becomes bimodal at 
a force ratio in the range 1.01 < 7 < 1.34. This observa- 
tion is in good agreement with the theoretical predictions 
of Section UlI Cl where a bimodal distribution is predicted 
for 7 > 7* = 1.257. 

To account for the dramatic changes evident in the 
measured distributions we have calculated p(x) for a 
periodically-driven oscillator. The exact shape of this 
distribution depends solely upon the force ratio 7, as 
noted in Section IIII CI Fixing 7 at the experimental 
value, we have evaluated numerically the series of Besscl 
functions given in Eq. 1231 The infinite series was trun- 
cated when convergence to four significant figures was 
achieved. The results of these calculations are plotted as 



the solid lines in Figure [TOl The agreement is remarkably 
good and is particularly encouraging given that there are 
no adjustable parameters in the comparison. 



C. Sensitivity Limit 

The smallest charge that can be detected by our tech- 
nique is fixed by the point where the peak in the spec- 
tral density at uj p (the arrowed peak in Fig. [8]) be- 
comes indistinguishable from the Brownian background. 
At this point the spectral power in a band of width 
Aoj centred on cu p is P S i g = Aoj ■ (fcgT '/ V^) / \lq 2 c + aj 2 ). 
The corresponding ratio of electric to Brownian forces 
is 7 2 = Aw • uj c /-k(uj1 + uJi) and the smallest detectable 



particle charge is therefore, 



e\Z?i l 



\Z12kBTrjaAuj 
E 



(26) 



where we have used Eq. 1201 In a discrete Fourier trans- 
form the minimum width Aoj of a frequency channel is 
determined by the total time t of data acquisition and 
equals, Aw = 2ir/to. In our system where a = 610 
nm, r\ = 1.38 mPa s, and the maximum field strength 
is E ~ 80 kV m _1 , the minimum detectable charge 
is about 0.25 e for an experimental run of to = 26 s. 
This is equivalent to a limiting electrophoretic mobility 
of - 2.5 x 10~ 12 m 2 V" 1 s _1 . 



D. Measurement of Charge Distributions 

We have shown that optical microelectrophoresis pro- 
vides an accurate, rapid and highly sensitive determina- 
tion of the effective charge on an individual trapped parti- 
cle. In many cases, of course, the population from which 
the single particle is selected is heterogeneous. There 
is a distribution of particle charges which reflects, for 
instance, the stochastic nature of the charging mecha- 
nism or variations in the surface chemistry or the size 
of individual particles. Many of the conventional meth- 
ods used to probe the electrokinetics of colloidal suspen- 
sions (microelectrophoresis, laser Doppler electrophore- 
sis, and electroacoustics) only allow for the calculation 
of the average mobility. They provide either little or no 
information on the distribution of particle charges. In 
many cases, however, knowledge of the charge distribu- 
tion is of vital importance. Colloidal suspensions have, 
for example, proved to be valuable model systems in con- 
densed matter physics. There is strong evidence that the 
width of the charge distribution (the charge polydisper- 
sity) significantly influences the glass transition and at 
high levels may suppress the freezing transition in these 
systems [2l|. It is therefore important that techniques 
are available which allow an accurate characterization of 
the charge distribution. 

We consider a suspension of colloidal particles with (ef- 
fective) charges Ze which are distributed according to a 
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normalized charge distribution function P(Z). The mean 
charge and charge polydispersity oz are defined by 



P(Z)ZdZ 



P(Z)(Z-Z) 



1/2 



(27) 
(28) 



To characterize the charge distribution we trap and mea- 
sure the effective charge of 50-100 individual particles. 
Although a little tedious these repeated measurements 
arc still quick. For instance it takes a couple of min- 
utes per particle to manually measure its charge: 20-30 
s scanning time to locate a particle, a short delay of 30 
s to center the particle trap and a further 30 s for data 
collection. 

To explore the range of applicability and resolution 
of the technique we determined the charge distribu- 
tion in nonpolar suspensions of PHSA-coated PMMA 
spheres. These colloidal spheres have been widely 
used in fundamental studies of nucleation, crystalliza- 
tion and glass formation since the coating with PHSA is 
thought to provide a hard-sphere type interaction. Fig- 
ure II If a) shows the charge distribution measured for 
PHSA-coated PMMA spheres suspended in clean dry do- 
decane. Clearly the particles are very weakly charged, 
with an average colloidal charge of just three electrons 
per sphere, equivalent to an electrophoretic mobility of 
fx = —3 x 10~ n m 2 s" 1 V" 1 and a zeta potential of -3.5 
mV. Indeed the particle mobility was sufficiently small 
to be undetectable on a commercial PALS instrument. 
The low mobilities however presented no problems for the 
optical microelectrophoresis measurements, as the data 
in Figure [TlTa) confirms. Although the charge on the 
PHSA-coated spheres is small, it is clearly non-zero and 
so we expect a weak, soft repulsion between the particles. 
To gauge the level of this repulsion we estimate the aver- 
age contact value Uq of the interaction potential between 
spheres. An average charge of Z = —2.9 equates to a 
value of U /k B T = 0.2 in dodecane (Eq. [J). Although 
the value is small in comparison to thermal energies it 
can not be ignored in quantitative studies. This is high- 
lighted by, for instance, the recent studies of Auer et al. 
[22| who have explored the effect of a very similarly-sized 
weak repulsion on the crystallization rate. They found 
that soft spheres crystallize some two orders of magnitude 
faster than a comparable system of pure hard spheres. 
Clearly, for quantitative studies, it is important to fully 
characterize the interaction potential and optical micro- 
electrophoresis offers a promising alternative to existing 
techniques. 

The data in Figure [TTT a) reveals a second distinctive 
feature. The suspension is amphoteric, in that it contains 
both positive and negatively-charged spheres. Close to 
two-thirds of the particles are negative and the remainder 
positive. While the charge distribution seems to be ap- 
proximately Gaussian in shape, the width is broad with 
a charge polydispersity of oz = 1.5 ± 0.3, which is 



surprisingly large. 

As a further example of the applicability of the opti- 
cal electrophoresis method, Figure [TTT b) and (c) show 
the effect of adding the surfactants, Zr(Oct)2 [b] and Na- 
AOT [c], at concentrations of 2 mM to dilute PMMA 
suspensions in dodecane. These surfactants form reverse 
micelles in nonpolar solvents such as dodecane. A finite 
fraction of these micelles spontaneously ionize at room 
temperature (23|. Particles then acquire charge either 
by preferentially adsorbing the micelles of one sign or by 
the dissociation of a surface group, with the subsequent 
ionic species being solubilized within a reverse micelle. 
While the micelle-mediated chargingof colloids has been 
known for at least the last 60 years [J] , the exact mecha- 
nism of charging in nonpolar media remains problematic 
and is not well understood. Figure [IT] reveals that ad- 
dition of Zr(Oct)2 causes the PMMA spheres to become 
strongly positively charged, while Na-AOT generates a 
large negative particle charge. Using the Hiickcl equa- 
tion [HI we calculate the equivalent zeta potentials to be 
+63 mV for the Zr(Oct)2 treated spheres and -60 mV for 
the AOT containing suspensions. A comparable value 
for the zeta potential of around -38 mV was found by Ki- 
hara et al. [25| for PMMA particles in cyclohexane with 
added AOT. While Hsu et al. observed significantly 
higher zeta potentials, of around -140 mV for PMMA in 
dodecane. In the case of Zr(Oct)2, data is more limited, 
although Croucher et al. [26| reports that poly(vinyl ac- 
etate) particles are charged positive by the zirconyl salt. 
Surprisingly, we find that the charge distribution of the 
zirconyl and AOT suspensions are very sharply peaked 
with charge polydispcrsities of 5.0 ± 0.5% for the AOT 
system and 5.8 ± 0.6% for ZrOEH. In both cases the 
charge polydispersities are comparable, allowing for the 
sampling errors, to the size polydispersity of the colloid 
(4.6 ± 1%). Detailed information on the charging mech- 
anism in these systems will be published elsewhere. 



V. CONCLUSIONS 

In summary, we have demonstrated a novel technique 
for the direct measurement of small colloid charges com- 
monly found in nonpolar suspensions. The charge is 
measured by confining an individual particle to a har- 
monic optical potential and by following the excitation 
of the particle by a sinusoidal electric field. The trapped 
particle forms a strongly damped oscillator. By using a 
weak optical trap and an intcrferomctric position detec- 
tor we have shown that single particle optical microelec- 
trophoresis (SPOM) is capable of very high sensitivity. 
Surface charges on the level of a few elementary charges 
can be detected on individual colloidal particles, with an 
uncertainty of about 0.25 e. Other techniques of measur- 
ing the charge on a particle are generally less sensitive. 
Our measurements are rapid and reproducible, with data 
on a single particle being recorded in approximately 30 s. 
When applied to multiple particles our technique yields 
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information on the distribution of particle charges rather 
than simply recording the average charge, as provided by 
most conventional electrokinctic techniques. These char- 
acteristics makes SPOM ideal for the characterization of 
nonpolar suspensions and may prove useful to study fun- 
damental phenomena in colloid chemistry, such as the 
relaxation of the double layer. For example, one might 
be able to detect discrete changes in the charge on a par- 
ticle, due to the time-dependent relaxation of the double 
layer. Work exploring these possibilities is currently un- 
derway. 
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FIG. 1: Simplified experimental set-up of optical microelectrophoresis apparatus. WP, A/2-wave plate; M, gimbal mirror; BS, 
beam-splitting cube; L, lens; ND, neutral density filter; O, oil-immersion objective; C, high NA-condenser; QPD, quadrant 
photodetector; FG, function generator. 




28 Jim 



(b) 




FIG. 2: (a) Experimental geometry for the optical microelectrophoresis apparatus (not to scale), (b) Optical cell. 
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FIG. 3: The dimensionless ratio 7 of electric and Brownian forces. The trapped particle has a charge eZ e g and is placed in an 
oscillatory electric field of amplitude E and frequency uj p . The Brownian forces (shown on the right-hand axis) are calculated 
assuming u> p <C uj c and a force constant of kn = 5 fN nm _1 . The shaded region indicates the noise floor where the spectral 
density of the electric and Brownian forces are equal at ui p . The boundary is calculated assuming a = 500 nm, 77 = 1.38 mPa 
s and to = 26 s. In the shaded region, the charge on the particle is too low to be detectable by the current technique. The 
limiting charge sensitivity is ~ 2 e and 0.2 e for fields of 10 4 and 10 5 Vm -1 , respectively. 
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FIG. 4: The balance between thermal and electric forces. The full lines show the dependence of the asymptotic probability 
distribution, p(x), on the ratio 7 of the electric and thermal forces. The inset plot shows the curvature k of p(x), evaluated at 
the origin x — 0. The probability distribution is double peaked for 7 > 1.257 where n changes sign. 
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FIG. 5: Time trace of the x- and y-coordinates of a charged particle in a sinusoidal electric field. The field was aligned along 
the x-axis and had a frequency of 67 Hz and an amplitude of 13.3 kV m _1 . Periodic oscillations are not visible in the measured 
trajectory, which is dominated by random Brownian fluctuations. The relative size of the periodic oscillations expected are 
illustrated by the lower curve in (a). 




FIG. 6: The measured autocorrelation function C(t). The sample was a PMMA particle in dodecane (containing 0.035 wt% 
free PHSA-g-PMMA copolymer). Data was recorded for 26 seconds at each field. The driving frequency ui p was 420.8 rad s _1 . 
Fig. (a) shows that at short delay times the measured autocorrelations collapse to a single exponential decay. The dashed lines 
in (a) are least-squares fits to Eg. 1141 with a corner frequency of 295 ± 7 rad s _1 . Fig. (b) reveals the periodic correlations 
found at long delay times. The curves have been vertically displaced for clarity. 
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FIG. 7: The field dependence of the force ratio 7. 7 was obtained by fitting the autocorrelations displayed in Fig. [6]to Eg. I14I 
The linear variation of 7 with the field E confirms that the charge on the particle is independent of the strength of the applied 
field. The solid line yields an effective charge of \Z c g\ — 14.5 ± 0.2. 
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FIG. 8: The field dependence of the spectral density The sample was a 610 nm PMMA particle suspended in dodecane 

with 0.035 wt% added PHSA-g-PMMA copolymer. The variance of the power spectrum was reduced by averaging together 64 
spectra, calculated from consecutive trajectories each containing 2 14 data points (to = 1-64 s) recorded on the same particle. 




FIG. 9: The frequency dependence of the force ratio 7. The amplitude of the electric field was fixed at E = 75.8 kV nT 1 and 
the frequency uj p varied. The solid line is a least-squares fit of log against log(l + (lo p /uj c ) 2 ) with a gradient of -0.98 ± 0.02. 
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FIG. 10: The probability distributions characterising the restricted diffusion of an optically-trapped PMMA particle in a 
solution of 100 mM Na-AOT in dodecane. The corner frequency of the optical trap was 25.0 Hz. The points give the measured 
distribution in electric fields of amplitude (a) kV m -1 , (b) 22.8 kV m -1 , (c) 45.5 kV m~\ and (d) 60.7 kV m" 1 . The sinusoidal 
electric field had a frequency of 17.5 Hz. The solid lines were calculated from the theory for a driven harmonically-bound particle 
(Eq. 1221) using experimentally-measured parameters. The comparison has no adjustable parameters. 
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FIG. 11: The charge probability distribution measured for PMMA particles in dry dodecane with (a) no added surfactant, (b) 
2 mM zirconyl 2-ethyl hexanoate, and (c) 2 mM sodium-AOT added. The full lines are Gaussian fits. 



